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Abstract 

Quantum interference between one- and two-photon absorption pathways allows coherent control 
of interband transitions in unbiased bulk semiconductors; carrier population, carrier spin polar- 
ization, photocurrent injection, and spin current injection may all be controlled. We extend the 
theory of these processes to include the electron-hole interaction. Our focus is on photon energies 
that excite carriers above the band edge, but close enough to it so that transition amplitudes 
based on low order expansions in k are applicable; both allowed-allowed and allowed-forbidden 
two-photon transition amplitudes are included. Analytic solutions are obtained using the effective 
mass theory of Wannier excitons; degenerate bands are accounted for, but envelope-hole coupling 
is neglected. We find a Coulomb enhancement of two-color coherent control process, and relate it 
to the Coulomb enhancements of one- and two-photon absorption. In addition, we find a frequency 
dependent phase shift in the dependence of photocurrent and spin current on the optical phases. 
The phase shift decreases monotonically from 7r/2 at the band edge to over an energy range 
governed by the exciton binding energy. It is the difference between the partial wave phase shifts 
of the electron-hole envelope function reached by one- and two-photon pathways. 
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I. INTRODUCTION 



The phenomenon of quantum interference can be used to control physical and chemi- 
cal processes. One approach, the 'ra + m' scheme, uses a two-color light field to interfere 
n- and m-photon transitions.ii*2i^ Interference between one- and two-photon transitions, for 
example, allows controllable polar asymmetry of photoelectrons in atomic ionization, ^'^ con- 
trollable dissociation of HD+,^ and controllable photocurrent injection in unbiased solids due 
to free carrier absorption,— impurity-band transitions,— and quantum well bound-continuum 
intersubband transitions.— In biased asymmetric semiconductor double wells, '1-1-2' inter- 
ference allows control of carrier population and THz emission.— Interband '1 -|- 2' interfer- 
ence in unbiased semiconductors, which is our interest here, allows independent control of 
electrical current injectio n^^i^^ and spin current injection . ^^i-^^i^^i^^i^^i^? Furthermore, in non- 
centrosymmetric semiconductors, it allows independent control of carrier populations (i.e., 
absorption) ^^1^° and carrier spin polarization . In each scenario, the experimenter can 
control the interference by adjusting the phases of the two colors. 

The controllable optical phases are not the only source of phase between the transition 
amplitudes. In general, there is also a material-dependent intrinsic phase.— Phenomenolog- 
ically, the intrinsic phase appears as a phase shift in the dependence of the process on a 
relative phase parameter of the optical fields. Additionally, selectivity between two processes 
is possible when their intrinsic phases differ^^; for example, '1-1-3' experiments on diatomic 
molecules have controlled the branching ratio of ionization and dissociation channels.— The 
intrinsic phase can be strongly frequency dependent near resonances,— and the hope that 
it might serve as a new spectroscopic observable^^*^ has led to efforts to understand its 
physical origin. 

Whereas a resonance is necessary for a phase shift to a '1-1-3' process,— it is not necessary 
for a phase shift to a '1-1-2' process. For example, an intrinsic phase in the '1-1-2' photoion- 
ization of atoms is predicted from the simple model of a delta function potential . ^^'^^^ 

Nevertheless, microscopic theories for the interband '1 -|-2' processes in bulk semiconduc- 
tors have until now predicted trivial intrinsic phases.— The photocurrent, for exam- 
ple, was predicted to be proportional to the sine of the relative phase parameter for all final 
energies.— Each of these theories use the independent particle approximation, in which the 
Coulomb attraction between the injected electron and hole is neglected. That approxima- 
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tion is expected to be good for final energies well above the band gap, since in this case 
the electron and hole travel quickly away from each other. However, close to the band gap, 
one generally expects to see deviations from the independent particle approximation. In 
the one-photon absorption spectrum, for instance, it is well known that the electron-hole 
attraction is responsible for exciton lines below the band gap, and for an enhancement of 
the absorption above the gap known as Sommerfeld or Coulomb enhancement^ 

The effect of the electron-hole attraction on one-photon absorption has been studied 
with various degrees of sophistication. On the one hand, modern ab initio calculations 
that include Coulomb effects have recently given very good quantitative agreement with 
experimental spectra ,^'^i^^i^^i^'^i^^ although at the cost of significant computational overhead. 
On the other hand, simple models of Wannier excitons can describe Coulomb effects near 
the band edge of many direct gap semiconductors. These excitonic effects have long been 
understood quahtatively on the basis of a simple two-band model in the effective mass 
approximation,— which is even quantitatively accurate for typical semiconductors.— 

Excitonic effects on nonlinear optical properties of bulk semiconductors have also been 
studied in the effective mass, Wannier exciton approximatio n^'''i^^i^^i^°i^'^i^^i^^i^"^i^^i^^i^'''i^^i^^i^° 
and only recently with ab initio methods. The two-photon absorption spectrum shows a 
different set of exciton lines and a Coulomb enhancement that is weaker than its one-photon 
counterpart. 

One- and two-photon absorption spectra have been measured sufficiently often that ex- 
citonic effects on them are well established. In contrast, semiconductor '1+2' interference 
experiments have been done typically at only a single energy and typically many exciton 
binding energies above the band-gap. Moreover, these initial experiments lacked an abso- 
lute calibration of the relative optical phase, and thus were insensitive to the intrinsic phase. 
Such a calibration is possible,— and could be used to verify the predictions we present here. 
A nontrivial intrinsic phase would have implications for the use of '1 + 2' current injection as 
an absolute measurement of the carrier-envelope phase of an ultrashort optical pulse .■^^'^^'^ 

In the present work, we extend the theory of '1+2' coherent control of bulk semiconduc- 
tor interband transitions beyond the independent particle approximation, employing a set of 
approximations that are valid close to the F point of a direct gap bulk semiconductor. Our 
investigation is limited to a perturbative treatment in the fields. In this limit of low pho- 
toinjected carrier density, the only inter-particle interaction of importance is that between 
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a single electron and hole. We show that, due to the electron-hole attraction, a nontrivial 
phase shift does in fact occur in the control of current and spin current, but not in the 
control of carrier population or spin polarization. The intrinsic phase can be understood 
in terms of partial wave phase shifts due to the Coulomb potential between electron and 
hole. In addition, we find an enhancement of each process, and relate it to the Coulomb 
enhancements of one- and two-photon absorption. 

In the next section, we establish notation necessary to describe the '1-1-2' processes in 
terms of one- and two-photon transition amplitudes. In section IIIIl we present the micro- 
scopic model. The transition amplitudes are worked out in section HVl The final expressions 
for the '1-1-2' effects are given in section |3 and numerical results for GaAs are presented. 
In section IVIl further understanding of the enhancement and intrinsic phase is discussed, 
and we examine the ratios often used as figures of merit for '1-1-2' effects. Intermediate state 
Coulomb enhancement is examined in Appendix ^ For reference, the current injection 
tensor is worked out for parabolic bands in Appendix |B1 



II. PRELIMINARIES 



The rate of photocurrent injection into an unbiased bulk semiconductor by a two color 
light field E(t) = Fi^exp{—iut) + E2uj exp{—i2u!t) + c.c. can be written 

f =^e:^e:^ei+c.c., (1) 

where superscripts denote Cartesian components and repeated indices are to be summed 
over.— The fourth rank tensor 7/(7), called the current injection tensor, describes the material 
response. It is purely imaginary in the independent particle approximation,— but can be 
complex in general. We define the intrinsic phase 5^^^^ of the component rj^^^^ as 

5''^' = arctan (rc (r^j^f /Im (^r]f^'^ ) (2) 

so that it is zero or vr in the independent particle approximation. When the electron-hole 
interaction is included in the set of approximations used below, all the components of 77(7) 
have the same phase. That is, 

(3) 



ijkl 



ijkl 
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The intrinsic phase S appears as a phase shift in the dependence of the current injection 
on the phase of the optical fields. For co-linearly polarized fields (E^ = Ei^exp{i(p^)^ and 
Fi2w = -^20) exp(202a;)x), for example, the current injection is 

^ = 2ElE2^ \T]^Jf^\ sin (20^ - 02.. - S) x. 

We are ignoring scattering processes through which the current would relax to a steady state 
value under continuous illumination, or would decay to zero following pulsed excitation. The 
current injection discussed here can be used as a source term in hydrodynamic equations 
that treat the scattering phenomenologicallyii*^^^ or in microscopic transport equations.— 
Coulomb effects other than the excitonic effects we consider here play a role in scattering, 
especially at high densities of excited carriers. Such Coulomb effects are outside the scope 
of this paper. 

The current injection tensor can be written in terms of one- and two-photon transition 
amplitudes. We take as the initial state a clean, cold semiconductor. In a Fermi's golden 
rule calculation for the ballistic current, light produces transitions to final states \n) with 
velocity v„„, probability amplitude c„, and energy hjn above the ground state. The final 
states will be specified in detail in the next section; here the label n represents the set of 
quantum numbers for either an interacting or independent electron-hole pair. Thus, 




n 



where e is the electron charge (negative), is a normalization volume, and f2n is the 
amplitude for an z-photon transition. The Qn^ take the form 

f^i^^ = E2. ■ D« (5) 
= E.E. : , (6) 

where the vector and second rank tensor D^^'' depend only on properties of the material. 
We consider them in detail in section HV1 

In the last equation of (0)), the first term is (one-photon) one-color current injection, the 
second term describes two-photon one-color current injection, while the interference term 
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describes the '1+2' current. Because of their different dependencies on the electric field 
amplitudes, these three terms can in principle be separated experimentally. But the first 
two terms vanish for centrosymmetric crystals, and the first vanishes even for zincblende 
crystals; we neglect them here. Excitonic effects on the first term were studied by Shelest 
and Entin .^^i^" The third term survives in all materials. Comparing its expression with the 
phenomenological form (^, we have 

= I? E <n {D^'*f {D^'P {2. - u^) . (7) 

n 

Even if scattering from impurities and phonons is neglected, the injection current de- 
scribed by (jH) does not capture the full current density; there are also optical rectification 
and 'shift' contributions to the current.— The one-color varieties of these have been stud- 
ied in some detail ,^°i^^i^^i^'^ but the '1+2' varieties have not. However, the different time 
dependencies of the three current contributions allows for their separate examination ex- 
perimentally, at least in principle, and rough order-of-magnitude estimates indicate that 
typically the injection current will be the largest; it is the only contribution to the current 
we treat here. 

As in (0}, the carrier injection rate can be written as 



d I ^ d 

ft 

dt L3 ^ dt 

n 



UJr. 



where n is the number density of injected electron-hole pairs. The first two terms in (jH)) are 
the usual one- and two-photon absorption rates, which we denote by 77-20^ and h^^ respectively. 
In terms of one- and two-photon coefficients Ci) and Cif that depend only on the properties 



'(1) ^''^ H2) 

of the material, they can be written as h2uj = ^(i)-E'L-^2i!> and = Cij^^ ElElE^ From 



^ = |EW'^*)^W^)'^(2^--«) (9) 

n 

^ = |EW^l"W^)"^(2---n)- (10) 

n 

The third term in (jHJ, denoted n/, allows population control as discussed and observed by 
Eraser et alM^ It can be written in terms of a third rank tensor ^^-^^ 

nj = ^f^E*JE*JEl + c.c., (11) 



where 

C = |EW'^*)"W'^)'^(2---")- (12) 

n 

It is purely real in the independent particle approximation . 

Expressions such as and (jHl) can also be written for carrier spin polarization and spin 
current.— The interference terms of these describe '1+2' spin control^ and '1+2' spin current 
injection,— which can be written in terms of material response pseudotensors C^^ — and 
^tjk^im^ig 2;ggpgQ|;iYgiy_ ^j^g independent particle approximation, C^-'j^'' is purely imaginary,— 
while /i^^)'*" is purely real.— 

The phases of the material response tensors ?](/), ^(/), and are related to the 
phases of the one- and two-photon matrix elements, JDn^ and Dn^ . The one- and two-photon 
matrix elements also appear, respectively, in the one- and two-photon absorption coefficients, 
as can be seen from and There have been many theoretical investigations of one- 

and two-photon absorption near the direct gap of bulk semiconductors that include excitonic 
effects . ^^1^^ However, since one- and two-photon absorption are insensitive to the phases of 
the transition amplitudes, those calculations took no care to get the phases of the transition 
amplitudes correct. In the next two sections we find the transition amplitudes with the 
correct phases including excitonic effects. 

III. MODEL 

We first review the two-band effective mass model of Wannier excitons; the two bands are 
nondegenerate conduction and valence bands that are parabolic and isotropic with a direct 
gap E^^ at k = (the F point). -^^^ It has been used to study excitonic effects on one-photon 
absorption,i^ two-photon absorption,^^ and other nonlinear optical processes . '^^i^^'^^'"^^'"^^'"^^'^' ' ' 
We subsequently describe a generalization that accounts for degeneracy and multiple bands. 
It has been used for two-photon absorption,— and has been implied whenever two-band 
results have been applied to actual semiconductors. 

The total Hamiltonian of the system can be written in the form H = Hb + He + i^int(^)- 
Here, Hq = Hb + Hq is the field-free Hamiltonian made up of the single-particle part Hb and 
the part due to the Coulomb interaction between carriers He- Using the minimal coupling 
Hamiltonian, the optical perturbation takes the form H[^t{t) = — (e/c) A(t)-v+e^yl^/ (2mc^), 
where A{t) is the vector potential associated with the Maxwell electric field and v is the 
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velocity operator associated with Hq. In the long wavelength limit, the position dependence 
of Ait) is neglected, and thus the second term in Hi^tit) may be neglected, since it can 
simply be absorbed in an overall time-dependent phase of the full system ket and hence 
cannot cause any transitions between states of Hq. Many approximate approaches to band 
structure calculation, including most pseudopotentials, and the truncation to a finite number 
of bands, implicitly assume an underlying field-free Hamiltonian that is nonlocal; there is 
then a correction to the interaction Hamiltonian i^int(t) in the velocity gauge .^^i^^ However, 
we neglect such nonlocal corrections, as has been the practice in previous calculations of 
coherent control effects.— 

The initial state is the "vacuum" |0); it corresponds to completely filled valence bands and 
empty conduction bands. If the Coulomb interaction were neglected in a two-band model 
consisting of valence (v) and conduction (c) bands, the final states would be of the form 
^ckO^k |0), where the operator al^^. creates an electron in an eigenstate of Hb, a Bloch state 
|?T,, k) with band index n and wavevector k. The photon momentum has been neglected, 
consistent with the long wavelength approximation. The Coulomb interaction couples states 
at different k; thus including He the final states are of the form 

\cVK)^J2'^cvi^)aLa,M, (13) 
k 

where k labels the state; its physical meaning is given below. In the effective mass Wannier 
exciton approximation, the Fourier transform 

C(r) = E^-(k)e^''^ (14) 

k 

which is the wavefunction for the relative coordinate between electron and hole, is a hydro- 
genic wavefunction satisfying 



V'^P:^{r) - \/(r)C(r) = {Ec. (ac) - E^J ^(r), (15) 

where fi~^ = + is the reduced mass in terms of the (positive) conduction and 
valence band effective masses, and V (r) is the Coulomb potential, V{r) = e^/ (er), screened 
by the static dielectric constant g .^OiSe^o rjj^g state has energy 

E,,{K) = -—-rEl,. 



'•cv 
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We choose the states to be normahzed over the volume by (m, k|n, k) = Sn,mSk,k' 
and {cvK,\cvK,') = 6,-^^,'] as a result ipcvi'^) is unitless, having the normalization 

Our focus in this paper is on the unbound solutions to (jl5|) : bound exciton states lack 
relative velocity between the electron and hole, and hence do not contribute to the ballistic 
current or spin current. For a Fermi's golden rule calculation of the current or spin current, 
the unbound state must behave asymptotically like an outgoing plane wave in the relative 
coordinate between electron and hole; k is the wavevector of the outgoing plane wave. 
Specifically, we must use "ionization states" rather than scattering statesj^ as was done for 
atomic '1+2' ionization.— They are related by {i'cvi'^))ion ~ [('^cti'*('^))scatt]*-^^ Calculations 
of one- or two-photon absorption are insensitive to an error in this choice of boundary 
condition, but the present calculation is not, since it is sensitive to the relative phase of the 
transition amplitudes. 

The ionization state wavefunctions that solve (|T5|l can be expressed as a sum over partial 
waves, 

oo Y (l + / ,■ \ 

^r.(r) = e^^Y. loi^r (2^«:r)^ e— (^) / + 1 + 2/ + 2; 2tKr , 

(16) 

where = e^^/ {p-cv^'^) is the exciton Bohr radius, and k, and r mean and |r[''^. The P/ 
are Legendre polynomials, i Pi is a confluent hypergeometric function, and F is the Gamma 
function. 

Such a two-band model of Wannier excitons is useful for the description of many optical 
properties. However, near the band gap at the F point of a typical zincblende semiconductor 
there are, counting spin, eight bands: two each of conduction (c), heavy hole {hh), light hole 
(Ih) and split-off hole (so). Other bands, especially the next higher conduction bands, can 
also be important for some processes, especially for population and spin control. 

The existence of multiple bands and band degeneracy modifies the exciton Hamiltonian, 
i.e., the operator acting on ipcvi^) i^ i^ft side of ()15|). In the effective mass approxima- 
tion, using a basis of F point states, the kinetic part of the Wannier exciton Hamiltonian 
has a matrix structure.—"^ Even though this approximation neglects band warping, non- 
parabolicity, and inversion asymmetry, the Hamiltonian lacks analytic eigenstates.>i^ This is 
essentially due to the degeneracy of the hh and Ih bands at the F point. As a result of the 
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difference between rrihh and mih there is 'envelope-hole coupling, — which is a spin-orbit-like 
coupling between the orbital angular momentum of the exciton envelope function and the 
total angular momentum of the valence band F point Bloch functions.— Baldereschi and 
Lipari split the effective mass Hamiltonian into a sum of terms based on symmetry, and 
showed that in a spherical approximation the envelope-hole coupling could be treated as 
a perturbation to the diagonal part, which has analytic, hydrogenic eigenstates.— >^ In or- 
der to extract the main physics, while preserving the simplicity of the two-band model, we 
neglect envelope-hole coupling entirely. In this approximation, ()15p remains valid for each 
conduction-valence band pair, however one must use 'average' effective masses for degen- 
erate bands. Specifically, the effective mass of the valence bands hh, Ih, and so is m/7iL, 
where m is the free electron mass, and 71^ is one of the Luttinger parameters.™ The upper 
conduction bands have a different average effective mass. Note that '?/'^(r) is independent 
of c and v within the set of bands {c, hh, Ih, so}. The effect of envelope-hole coupling has 
been studied for exciton bound states,'^2i22i2S but not for optical processes involving unbound 
excitons in the continuum. 

Even within this model, the presence of multiple bands provides two types of terms in the 
sum over intermediate states in the two-photon amplitude: two-band terms, in which the 
intermediate and final states are in the same exciton series [i.e., two states of the form ()13|) 
with the same c and v] , and three-band terms, in which the intermediate and final states are 
in different series. Three-band terms are important for some processes but not for others. 
For current control, three-band terms are important for cross-linearly polarized fields,-^ and 
for spin current control, they are important for the spin current due to co-linearly polarized 
fields.— Three-band terms are essential for population and spin control . ^^1^° 

The velocity matrix elements involving the state \cvk,) are 

(ct;..|v|0) = 5^(AS,(k))*v,4k) (17) 

k 

{cvK,\ V \c'v'^') = J2 (Ka (k))* A^:^, (k) (v,,, (k) - V,,, (k) S,y) , (18) 

k 

where \nm (k) = {n, k| v |m, k) is the velocity matrix element between Bloch states. 
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IV. TRANSITION AMPLITUDES 



In the independent particle approximation, the transition amphtudes are 

^^i^-r^ = ^^E2. ■ V,. (a.) , (19) 



and fi^'i'"'^ = J2c',v' ^cc'vv'l where 

(^(2-free) ^ f 6 >^ ^Eg; " (Vcc' («^) - ^c,c'V^'^ («^))) (Eg; ■ Vg/^/ (k)) 

With excitonic effects included, using the perturbation iJint it) to second order gives the 
transition amplitudes 

fiill = ^E2.-(ct;A^|v|0), (21) 



and 

o(2) = (iS V (E^ • {cvn\ V \c'v'n')) (E^ ■ {c'v'k'\ v |0)) 



(22) 



where the sum over intermediate states is over both bound and free excitons. The two- 
photon transition amplitude is more difficult to deal with due to the sum over intermediate 
states; however, in our set of approximations it has been done exactly .^^i^^i^" 

In order to proceed analytically, it is common to use (fTTj) and ()18|) . and then make 
an expansion in k of the velocity matrix elements v„,m (k) about the T point 
However, due to the degeneracy at the F point, the coefficients of such an expansion can 
depend on the direction of k.^^ To proceed, we note that Wannier excitons have large spatial 
extent and hence only a small region of wavevectors is important for them, i.e., A'^^ (k) is 
peaked in the region of k near k.— This is especially true for final states with energies above 
the band gap. Thus, we expand \nm (k) about the F point, approached in the direction k, 

<m (k) = <^ (^) + k ■ Vk<„(/^) + . . . , (23) 

where (k) = limA->o F| v |m, Xk) and Vkf^^ (k) = liniA^o (?^, F| v |m, Xk). 

Optical transitions due to the first term in (j^Hj) are called 'allowed', while those due to the 
second term are called 'forbidden'. We restrict ourselves to materials for which the 'allowed' 
valence to conduction band transition does not vanish. Keeping only the 'allowed' term in 
(HZj),^ 

(cr;/.! V |0) = (/^) (C(r = 0))* . (24) 
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For the intravalence and intraconduction band transitions, the first two terms of ()23p in ()18|) 

d^7 



/d/ I I 
-jj ii^cv (r))* 

— (C.(r))* V^X'(r). 

In particular , ^^1^^ 

{cv^\ V \cvH,) = J2 (K. m* (k) — k = -I— [ ^ (C.(r))* VrC,(r) = 



(26) 

For Ge and for simple models of zincblende semiconductors that neglect lack of inversion 
symmetry, the first term in (j25j) always vanishes. This means that there are only 'allowed- 
forbidden' two-photon transitions (the interband transition is allowed, while the intraband 
transition is forbidden). When the first term is nonvanishing, there are 'allowed-allowed' 
two-photon transitions. In principle, for materials that lack a center of inversion, there is 
also a small contribution to the 'allowed-forbidden' two-photon transition from the first term 
in and the term in {cvn\ v |0) that comes from the second term in (j2Sl); we neglect it in 
what follows, but note that when compared to the dominant 'allowed-forbidden' contribution 
that we consider here, it has a different Coulomb enhancement but the same intrinsic phase 
(see Eq. 2.32 of Rustagi et al.)M 

We write fii^li = 0,''^'^^ and discuss the 'allowed-forbidden' and 'allowed-allowed' 

transitions separately. 

Using (fTI)|) . and (j2H), the one-photon transition amplitude isM 

where only the 'allowed' term is kept in Q^Juk^^\ The transition is to an s-wave. The 
one-photon absorption coefficient is proportional to the norm of ficiL [see dHj) and ©])•" 

For 'allowed-forbidden' two-photon transitions, substituting (j24j) and the second term of 
(123) into (1221), 



^^^"^ = E ■ (^)) i^^y^^^c i^) - 4,c'Vk<,„ (/^)) ■ M,,,w (k) , (28) 

c'v' 

where 



M,,,,,, (k) = -iJ d\ (V-r. (r))* VG,'.' (r, 0; - Ef,,,) , (29) 
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and we have used the Coulomb Green function, 

which is known analytically.'^^ In particular, 

G„(r. 0; fi. - E?„) = (1 - . (-il- ) , (30) 



where we define 



Ten — 



5, 



with the integral representation 



Bcv = / (2/ici;Oct,) is the exciton binding energy, and Wy,i/2 [z) is a Whittaker function 



Since the Green function depends only on the magnitude of r, only the p-wave of the final 
state survives the integral in Eq. ()29|) over the angles of r, j dQPi (f ■ k) f = Ank/S. The 
integral over r can be done using^— 



r'^-'e-P\F, («; a; \r) dr = T (a) / (32) 

(p- A) 



The final result is 



where only the 'allowed-forbidden' term is kept in ^f^'m''!^ ^^"^ 



ri S (jif )'^"'''' exp ( --^ arctan {ac'v'K.lc'v'S) 
N^l (-) - (1 + aL'-^7.\0 2 / -^^^^ r; , '-dS. (34) 

For 'allowed-allowed' two-photon transitions, substituting (j24j) and the first term of (|25p 
into (1221), 

= (£) ' E ■ ['^-'^^c (z^) - ^c,c<>. (^)] ■ V,,,, (/i) 
X /i y (V^,^, (r))* G,,,, (r, 0; - Ef,„,) . 

Since Gc'v' depends only on the magnitude of r [Eq. (jHOI)], only the s part of the final state 
will survive the integration over angles of r. Again we use (|HT|) for the Whittaker function. 
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The integral over the magnitude of r can be done using an identity obtained by taking a 
derivative with respect to p of both sides of (jH^ . Finally, 



"Sr' - (2^) - i) E "-l^Sivi?^. («) , (35) 



where only the 'allowed-allowed' term is kept in ^^^cWk' ^^"^ 
iVi%"), = (l + (a,v«:7cv)')2 



1 _l_ ^-x 'Tc'v' exp ^— arctan {ac'v'i'^Jc'v'S)^ 



X 5 fl - 5^^^^^') ( — ^ ) ^^^1: , '-dS. 

Jo \ acv J\1-SJ (1 + (a,v/^7cv^) ) 

(36) 

This agrees with Eq. 2.28 of Rustagi,!^ but note that we have defined N^^'vl' ~ 
[1 + {ac'v't^lc'v')'^) Is,k {1^), where Is,k{i^) is given, with a typographical error, in Eq. 2.25 
of that paper. 

The factors N^^'vv' ^cc'vl' ' which appear in (jH^ and (IHUj) are the enhancements due to 
the Coulomb interaction in the intermediate states; they are discussed further in Appendix 

El 



V. RESULTS 



The one- and two-photon transition amplitudes were presented in the previous section on 
the basis of an expansion in k of the Bloch state velocity matrix elements. The 'allowed' one- 
photon transition amplitude fi^^-* is in ()27p. the 'allowed-forbidden' two-photon transition 
amplitude is in (jH^ . and the 'allowed-allowed' two-photon transition amplitude is in (jSHJ. 
From them, DcJ,lt and DcS^ may be extracted by comparison with the definitions in (0) and 
©■ 



A. Current injection 

The '1-1-2' current injection is dominated by interference of 'allowed' one-photon tran- 
sitions and 'allowed-forbidden' two-photon transitions.— Substituting these into and 
using the Gamma function identities F (a; -|- 1) = xF (x) and 

F (1 — ix) F (1 -|- ix) = TTx/ sinh (ttx) , (37) 
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yields our final result for the current injection tensor 



4) = E f 1 + rV) " ^^-'^-^ ^ii' M ^L. (38) 

c,v ^ ' d ,v' 



where 



1 2hu-E? 



and 



V, 



^) = ■ W ( = — (1 - exp (-27r/a;)) , 40 
sinh(7r/xj x 

27re . , . / ^r^-fronuX-''^ 



?i' - ^ E (-) (^iewt)' (^^"i^^^^)' ^ (2-^ - ^™ //^) , (41) 
with (k) = (vcc (k) - v^^ (k)) = hnjii^^ and 

/p,(2-free)V'^ _ ^ \ ^ {(Vcc' (/^) (k)),Vcv(k)V 

where {vi, V2}*'^ = {v\v2 + v{v2)/'2. and 

[Dt-'-^)' = i^v[^{^). (43) 

Note that only the 'allowed' part of PHj) and the 'allowed-forbidden' part of (jl^ should be 
retained for a consistent solution. We have written (j38p to separate the parts due to the 
electron-hole interaction. In the independent particle approximation, the current injection 
tensor r/J^^J^^^) isii 

ijkl \ ^ ijkl (AA\ 

^(/-free) ~ 2-^ Vcc'vv''^ \^^) 



c,c'v,v' 



it is evaluated for parabolic bands in Appendix ^ 

For GaAs, we present in Fig. [T]the magnitude of rj^j-j^^, based on r]^§^^, calculated by two 
methods. The first method, described in Appendix El uses isotropic parabolic bands and 
includes only two-band terms; it uses effective mass ratios for conduction, heavy hole, light 
hole, and split-off bands of 0.067, 0.51, 0.082, and 0.154 respectively, Ep = 27.86 eV, the 
fundamental band gap Eg is 1.519 eV, and valence band spin-orbit splitting is 0.341 eV .^^i^^ 
The second method solves the 8 x 8 k ■ p Hamiltonian including remote band effects, but 
in a spherical approximation with warping and spin-splitting neglected by replacing 72 and 
73 with 7 = (272 + 373) /5^^; the calculation is nonperturbative in k (hence it includes 
band nonparabolicity) and it includes both two- and three-band terms in the two-photon 
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Excess photon energy 2hu! — Eg (meV) 

FIG. 1: Magnitude of the diagonal element of the current injection tensor for GaAs with [Eq. ()38() ] 
and without [Eq. (|44() ] excitonic effects. The grey dotted and dash-dotted lines are based on a 
parabolic band calculation of rj^^^^i that only includes two-band terms; the dotted line includes 
excitonic effects, while the dash-dotted line does not. The black solid and dashed lines are based 
on r]^^,^^i calculated with a nonperturbative solution of the 8 x 8 k • p Hamiltonian; the solid line 
includes excitonic effects, while the dashed line does not. 

amplitude.~ The solid and dotted lines in Fig. [T] are calculated with (|HHj). and hence include 
excitonic effects; the Coulomb enhancement part of the calculation uses B^v = 4.2meV— 
and the band parameters listed above. Note that the solid black line in Fig. ^ is inconsistent 
in the sense that the Coulomb enhancement is based on an expansion in k, whereas the 
free-particle result that it enhances is nonperturbative in k; nevertheless, such an approach 
has given good agreement with experiments for one- and two-photon absorption . ^^'^7 

The Coulomb enhancement of //(/) can clearly be seen in Fig. ^ There is a kink in each 
curve at excess photon energy 341 meV corresponding to the onset of transitions from the 
so band. At higher energies, the Coulomb enhancement of so transitions is larger than 
the Coulomb enhancements of hh and Ih transitions, since the former transitions are to 
conduction band states with lower energy. Hence, the kink in 77(7) is enhanced by excitonic 
effects. 

We extract the intrinsic phase of rj^^^'^ using Q . The solid line in Fig. |21 is the intrinsic 
phase of r]^^^^ calculated for GaAs with the nonperturbative 8 x 8 k • p Hamiltonian band 
model; the result for the parabolic band model is nearly identical. Since we have used a 
spherical exciton model, the intrinsic phase is the same for all components of ff^pl'- The 
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FIG. 2: Phase shift of the current (intrinsic phase of mi^f-^^) in GaAs due to excitonic effects. The 



sohd hne is calculated with Eqs. ((3) and and the dotted line is calculated with Eq. 1)45(1 . The 
inset is Eq. (|45() plotted in scaled units. 

intrinsic phase has its maximum value of 7^/2 at the band edge, and goes to zero as the 
light frequency increases. The decrease is smooth except for a small kink at the onset of 
transitions from the so band. In fact, for excess photon energies less than the split-off energy, 
the intrinsic phase has the simple analytic form 



Equation ()45|) is plotted as the dotted line in Fig.|21 compared to the solid line, it is identical 
below the onset of so transitions, and it makes a good approximation above the the onset of 
so transitions. Since only depends on the excess photon energy scaled by the exciton 
binding energy, we plot it as a function of this scaled energy in the inset of Fig. |21 it is useful 
for finding the intrinsic phase of materials other than GaAs. 

In 77(7) [Eq. (jHH|) ]. the two- and three-band terms have different intermediate state Coulomb 
enhancement N^^'vi'- -^^^ many materials, however, N^^^j^, is approximately equal for all the 
terms rf^^l^t that contribute significantly to the total ?7(/^/j,ee)' shown in Appendix 1X1 for 
GaAs. Thus, at photon energies for which transitions from the heavy- and light-hole bands 
dominate ?7, the Coulomb enhancement becomes approximately independent of the sum over 
bands and we can make the simplification 




(45) 




(46) 



where the intrinsic phase is given by and 




(47) 
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FIG. 3: Approximate Coulomb enhancement factors. The solid line, applicable to current and 
spin current control, is F^_j [Eq. (|47jl ] with Nccvv = and the dotted line, applicable to carrier 
population and spin control is F^l [Eq. ^] with Ntvv = 1 

The Coulomb enhancement factor F^^^ (uj) is plotted in Fig. |21 with the approximation that 
Nccvv = 1 (see Appendix 

B. Carrier population control 

The '1+2' carrier population control is dominated by interference of 'allowed' one-photon 
transitions and 'allowed- allowed' two-photon transitions . Substituting these into (|12j). 
and using the Gamma function identity (j37|) . we find 

^ = E ^ («-'^-) E M et', (48) 

c,v c'v' 

where 

ew = ^ E (^2wt (Diir-'y^ (2- - Ecv m (49) 

K 

and and Dct^'^'^^'* are given by and (jlHjl . Note that only the 'allowed' part of 

PHj) and the 'allowed-allowed' part of (jl^ should be retained for a consistent solution. In 
the independent particle approximation, 

^(Ifree) ^ E ^cl'ra' • (50) 

Thus, population control has a Coulomb enhancement due to excitonic effects, but no phase 
shift. 

Note that (pSj) gives the population control tensor at final energies above the band edge. 
There can also be population control of bound excitons when both one- and two-photon 
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transitions are to the same excitonic state. This can occur, for example, at s excitons due 
to allowed-allowed two-photon transitions^ interfering with allowed one-photon transitions. 

If N^^l, (kcv) is approximately the same for all the terms that significantly contribute 
to .^(7), then, at photon energies for which transitions from the heavy and light hole bands 
dominate the Coulomb enhancement becomes approximately independent of the sum 
over bands, and we can make the simplification, 

~ -^a-a'?(/-free)' V'^^) 

where 

Fii) (u) ^ E (a^^K,,,) Nti! {kc.) . (52) 

The Coulomb enhancement factor F^^l (uo) is plotted in Fig. El with the approximation that 
Nccvv — 1 (see Appendix EJ. 



C. Spin current injection and spin control 

The '1+2' spin current is dominated by interference of 'allowed' one-photon transi- 
tions and 'allowed-forbidden' two-photon transitions, whereas '1+2' spin control is domi- 
nated by interference of 'allowed' one-photon transitions and 'allowed-allowed' two-photon 
transitions.— Under the approximations that led to ()46|) and ()51|). the spin current injection 
pseudotensor is 

f^S'^ = F!^hwi^s)f/^SZr (53) 

where F^^^ is given by (jUj), S is given by (jlH)), and /i(/-free) is the spin current injection 
pseudotensor in the independent particle approximation. Under similar approximations, the 
spin control pseudotensor is 

where Fa-1 is given by (j^ . and Cc^-free) is the spin control pseudotensor in the independent 
particle approximation. Spin control, like carrier population control, has a Coulomb en- 
hancement but no phase shift. There can also be spin control of bound excitons, but it has 
not been included in (jH^ . 
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VI. DISCUSSION 



We now examine the relationship between the Coulomb enhancements of the '1+2' pro- 
cesses and of one- and two-photon absorption; the latter are denoted by F^^^ and F^'^^ so 
that for i G {1, 2}, n^*) = n^ll^F'^^K The relationship is particularly simple at photon energies 
for which transitions from the heavy- and light-hole bands are dominant and intermediate 
state Coulomb enhancement is the same for each significant term in the sum over interme- 
diate states. The Coulomb enhancements for the '1+2' processes are then given by ()47|) 
and ()52|) . For one-photon absorption, F^^^ = S (acvKcv)-— In noncentrosymmetric semicon- 
ductors, two-photon absorption is dominated by allowed-allowed transitions just above the 
band gap, and by allowed-forbidden transitions at higher final energies; the cross-over point 
in GaAs is a few meV above the band gap.^^ At photon energies for which allowed-allowed 
transitions dominate two-photon absorption, from 

= E (a,„«:,„) (Arj-^) , (55) 

and thus 

= VF(i)F{2) and F^!^ = CVFWFi^\ (56) 

where C = (^Nccvv (ncv) /Nccvv incv)J 1 + iacvUcv)^"^, while at photon energies for which 
allowed-forbidden transitions dominate two-photon absorption, from (j33|l . 

= E (ae„«e.) (1 + (ae./€c.)"') {Nt'^ (/€,„))' , (57) 

and thus 

= (1 /C) VfWF(^) and F^^^ = y^pm^. (58) 

Note that, based on Appendix El C* ~ a/1 + -B™/ {2flLU — Eg), which is the ratio of the 
two curves in Fig. 01 In centrosymmetric semiconductors, there are no allowed-allowed 
transitions, and only (j3H|) applies. 

The '1+2' processes are often described by ratios. For example, a useful quantity to 
describe the current is the swarm velocity,— defined as the average velocity per injected 
electron-hole pair 

_ (dJ/dt) 

Vswarm = 71 TTTV' 

e [dn/dt) 

The swarm velocity is a maximum when the relative intensities of the two colors are chosen 
such that h2uj = riuj] returning to (jlj, if one associates the one- and two-photon amplitudes 
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with the arms of an effective interferometer, this condition corresponds to balancing that 
interferometer. For fields co-linearly polarized along x, the maximum swarm speed is 



1 



C (XX (XXXX 



(59) 



^{1)^(2) 

A useful quantity to describe pure spin currents is the maximum spin separation distance^^; 



it is proportional to /i(7)/ a/^(i)^(2)- As a consequence of (|^. the maximum swarm speed, 
and the maximum spin separation distance, are unaffected by excitonic effects when allowed- 
forbidden transitions dominate two-photon absorption.ifiP However, close to the band edge, 
where allowed-allowed transitions dominate two-photon absorption, excitonic effects increase 
these ratios by a factor C over their value in the independent particle approximation. In 
contrast, as a consequence of (j58p . excitonic effects do not affect the maximum control ratio 
for population and spin control {^{i)/ a/^(i)^(2) and C(/)/ a/^(i)^(2) respectively close 
to the band edge and decrease them by a factor C at higher photon energies for which 
allowed-forbidden transitions dominate two-photon absorption. 

In the terminology of Seideman,— the excitonic phase shift of the '1-1-2' current and spin 
current is a direct phase shift. It is due to the complex nature of the final state as it appears 
in the transition amplitudes. Thus it can be understood in terms of the partial wave phase 
shifts of the final state caused by the Coulomb potential between electron and hole. The 
Coulomb interaction is rather unique due to its long range nature, so we first suppose the 
potential between the electron and hole falls off more rapidly than 1/C. In that simpler 
problem, the final state wave function is written as 

^ — ' r \ rn / 

1=0 

where the u^,/ (r) are real.— If the potential between the particles is ignored, then the 
partial wave phase shifts, 6i (k) are zero. The allowed one-photon pathway reaches an 
s wave, while the allowed-forbidden two-photon pathway reaches a p wave. Substituting 
this form for the wave function into the one- and two-photon transition amplitudes, yields 
= Qk'^"^^^^ e^^"^'^^ fo (k) for the one-photon rate, where /o (k) is real and depends on 
Uk,o and 0,^'^'^ = ^n'^^'^^^e'^^'^'''^^ fi (k) for the two-photon rate, where /i (k) is real and 
depends on m^j q (f) and m^,! (r). Here ^1^'^^'^'^^ is the z-photon transition amplitude when the 
potential between the particles is ignored. It is then straightforward to see from ^ that the 
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relative shift of the partial waves is responsible for the phase shift of the current and spin 
current. That is, 

6 = 60-61. (60) 

The use of ionization states as opposed to scattering states was important to get the correct 
sign of the intrinsic phase. With scattering states, one would find 6 = 61 — 60. 

Due to the long-range nature of the Coulomb potential, the partial wave phase shifts 
have a logarithmic r dependent part, but it is the same for all partial waves and thus does 
not appear in the relative phase. The part of the Coulomb partial wave phase shift 61 (k) 
that does not depend on r is arg (F (/ + 1 + ^/(ac^,/^)))''^■^^; when inserted into (|6Up . this 
reproduces (j45j) . In contrast, the allowed-allowed two-photon pathway reaches an s wave 
and thus there is no phase shift for population control or spin control. 

The expression ()60p for the intrinsic phase in terms of the scattering phases is particularly 
simple, since each pathway connects to only a single parity. This contrasts with '1-1-2' 
ionization from an atomic s state, for which the one-photon transition is to a p wave and 
the two-photon transition is to both s and d waves; the intrinsic phase is thus a weighting of 
the p-s and p-d partial wave shifts.-2fi Materials for which the first term in (|23|) is forbidden 
(CU2O is an example) have these same selection rulesi^>^Si2£; hence, they will have an intrinsic 
phase with a similar weighting. 

The absence of a phase shift in population control can be connected to a symmetry of 
the second order nonlinear optical susceptibility. From considerations of energy transfer and 
macroscopic electrodynamics, ^(7) is related to the nonlinear susceptibility x^'^^ by 

= {leo/h) [x^'^'^^ {2u;- -u, -u) - x^^^'""' (-^; 2o;, -uj)] . (61) 

In the independent particle approximation,^^ 

X'-'^'^' {2u- -uj, -u) = [x^'^^^' {-uj; 2uj, -uj)] * , (62) 

which is a generalization of overall permutation symmetry to resonant absorption. As a 
result of (jnSI), Fraser et al. showed that ^(/) is proportional to Imx'^^^ and is thus purely 
real.i^ Our result that ^(/) remains real when excitonic effects are included suggests that 
()62|) holds more generally. In fact, it can be shown that (jU^ holds for any Hamiltonian 
symmetric under time-reversal so long as hw is not resonant. 
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VII. SUMMARY AND OUTLOOK 



We have extended the theory of interband '1+2' processes in bulk semiconductors to 
include the electron-hole interaction. Following previous theories, ^ii^i^'i^'i^iSi we have used 
a framework based on (i) a separation of the initial carrier photoinjection and the subse- 
quent carrier scattering, and (ii) a perturbative expansion in the optical field amplitudes, 
with injection rates obtained in a Fermi's golden rule limit for the bichromatic field. The 
injection rates for carrier population control, spin control, current injection, and spin cur- 
rent injection, have been described phenomenologically by tensors ^(/), V(i)^ /^{/)) 
respectively.— Like previous theories, we have used the long- wavelength limit, and 
neglect nonlocal corrections to the interaction Hamiltonian. But whereas previous theories 
of '1+2' photoinjection used the independent particle approximation, we have included ex- 
citonic effects. We have shown that excitonic effects cause (i) an enhancement of each '1+2' 
process, and (ii) a phase shift for current injection and spin current injection. Our main 
results, the modifications of the aforementioned tensors relative to the independent particle 
approximation are given in (jiUI) . (jHH), (j^ . and (jSH). These particularly simple results 
are valid at photon energies for which transitions from the heavy- and light-hole bands are 
dominant; more general results are given for 77(7) and ^(7) in (|nHj) and PHj). 

Our results are based on the effective mass model of Wannier excitons; degenerate bands 
are included, but we use a spherical approximation to the exciton Hamiltonian, and we 
neglect envelope-hole coupling. This is a good approximation for many typical semiconduc- 
tors, including GaAs, since the electron-hole envelope function extends over many unit cells 
due to the screening of the Coulomb interaction by the static dielectric const ant .i^2i22iZL22iZ^ 
As a consequence of making the spherical approximation, the phase shifts and Coulomb 
enhancements we find in this paper are independent of crystal orientation. 

Also, our results are limited to low excess photon energy since (i) the Wannier exciton 
Hamiltonian assumes parabolic Bloch bands, and (ii) we have truncated the expansion in k 
of the Bloch state velocity matrix elements, which is the basis of the transition amplitude 
expansion. By comparing the black dashed line and grey dash-dotted line in Fig. ^ one 
sees that higher order terms in k (for both bands and velocities) are important in GaAs for 
excess photon energies greater than about 200 meV. This can then be considered the limit 
of validity of our calculation. However, combining the Coulomb enhancement calculated 
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assuming parabolic bands with the nonperturbative independent particle approximation 
result (as was done for the solid black line in Fig. ^ likely gives a good approximation for 
a few hundred more meV; this was the case for one- and two-photon absorption .^^i^^ 

It is interesting to ask if there are other sources of intrinsic phases to the current (or spin 
current) besides the one that we have identified here, as these may produce spectral features 
in the intrinsic phase. One possibility is the coupling between bound so-c excitons and the 
unbound hh-c or Ih-c excitons, since it is known that the intrinsic phase can show spectral 
features near a resonance.- Another possibility is the envelope-hole coupling between the 
continua of unbound hh-c and Ih-c excitons that was neglected in our treatment. The effect 
of coupled continua in the general theory of 'n + m' phase shifts has not been considered to 
date. 

Finally, we note that the intrinsic phase and Coulomb enhancement may be greater in 
reduced dimensional systems, which have greater exciton binding energies. The carrier- 
carrier Coulomb interaction was included in the theory for '1 + 2' control of electrons in 
biased asymmetric quantum wells, although the intrinsic phase was not studied.-^ 

APPENDIX A: INTERMEDIATE STATE COULOMB ENHANCEMENT 

Consider the functions A^(^"f) and N^^~'^\ which appear in (j34p and (j36|) : we refer to them 
collectively as A^. First, note that due to the energy conserving delta function in ((Zj), k 
will be equal to k^v [see Eq. and thus is a function only of u, E^^, E^,^,, fi^v, and 

Hc'v'- Second, note that A^ is defined so that if the electron- hole attraction is turned off, for 
example by letting e oo, then A^ 1.^- This allows A^ to be identified as part of the 
Coulomb enhancement. In particular, A^ is the enhancement due to the Coulomb interaction 
in the intermediate states; if the Coulomb interaction is neglected for the intermediate states, 
A^ = 1.— [Note that Lee and Fan^ did not allow for v' ^ v in N (related to Jj in their 
notation).] 

Since the integrand is smooth for the parameter range of interest, numerical integration of 
A^ is straightforward; however, it need not be undertaken. Further simplification is possible 
since the parameter 7 can be considered to be much less than one. Since most materials 
have an exciton binding energy that is much smaller than the band gap, hw is detuned 
from the band edge by many exciton binding energies at photon energies consistent with the 
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Excess photon energy 2hL0 — Eg (meV) 

FIG. 4: The factor A^^*^"^) [intermediate state Coulomb enhancement, see (|33() ] for GaAs. The first 
(second) subscript is v {v'); subscripts for c = c' = Tq are not shown; Tg denotes the heavy and 
hght hole bands, and denotes the split off band. 

approximations made here. In GaAs, for example, when 2hw is within 500 meV of the gap, 
7 is at most 0.09. An expansion of A^^'^"^) for small 7, 

N^'-'^ = 1 + lie'.' + (^ln2 - 7,%, + (^So - ^cil^^ 7L' + O {it'.') , 

where 5*0 ~ 0.5633, shows that N^^~^ is approximately 1 and nearly constant as a function 
of to. The same is true of N^^~^\ which has the expansion 

Ar(-a) = 1 _ A (a^,,, - a,„) P + [i^^,) , 

where, with Si ~ 1.645, 

P = 7,,„, + 21n2 7^,„, + — i^a^'^'i^ + Si iX'- 

In fact, when = fJ'c'v', N^'^'^^ = 1 even to fourth order in ic'v'- Fig. HI shows a numerical 
integration of N^'^'^^ using the parameters of GaAs. 

APPENDIX B: EVALUATION OF THE CURRENT INJECTION TENSOR 

The tensor rf^^l^,, defined in fKTj) . can be used to calculate the current injection tensor 
with or without excitonic effects using (|38|) or (j44j) . It can be evaluated analytically in the 
approximation of parabolic bands. Part of the result for the 8 band Kane model has been 
given before but without the split-off band as an initial or intermediate state.— We here give 
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more detail, but only for the two-band terms. We denote the bands by a double index n and 
s, where n is one of {c, hh, Ih, so}, and s runs over the two spin states for each band. Since 
the Coulomb corrections to 77(7) in (jHS|) do not depend on the spin index, we can include 
the sum over spin indices from (j38|) in ff^^^i- And, since we are only calculating two-band 
terms, we set c' = c and v' = v. Thus 



Iccvv 



ire 



-v'5{2uj-E,,{k) /h) 



ire 



1 kl 



• ^ J , 



s' + (i ^ k) 



where 



2/^0 



{2hL0 - E, 



For V = Ih or hh, E^^ = Eg, while for v = so, E^^ = Eg + A, where Eg is the fundamental 
band gap, and A is the spin-orbit splitting. The interband velocity matrix elements are 
approximated by their value at A; = 0, but still depend on the direction of k. In terms of 
the orthogonal triple of unit vectors k, I and m, 

k = sin 6 cos </)x + sin 6 sin </)y + cos dz 
1 = cos 6 cos 0x + cos 6 sin 0y — sin 6z 
m = — sin 0x + cos (py. 



these matrix elements are: 



'c,s;hh,s' 



'c,s;lh,s' 



' c,s;so,s' 



2 V m 



1 Ep_ 
'2 V 3m 



- s,s' 

2kl + i\ay — irhax 



Ep 
6m 



kl — il(Ty + irhax 



where Ep is the Kane energy.iS^ Here, I is the 2x2 identity matrix and CTj are the 
Pauli spin matrices. Of course, for parabolic bands, the intraband matrix elements are 
(n, s, k| V \n, s', k) = 6s,s'^hk/mn, where m„ is the effective mass of band n. (In the proper 
Kane model, the effective masses are given in terms of the parameters Eg, A, and Ep, but 
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we treat them as additional parameters, which is equivalent to including remote band effects 
on the effective masses.) The sums over spin then yield 



^ hh,s' ■,c,s'^ c,s;hh,s' 



E 

s,s' 

E 

s,s' 



2m * ' 



Ep 
2m 



k^k' + -5k I 
3 ' 



Ep 

Ok I- 

3m ' 



The remaining angular integrals can be done using 



dnk'k^ = —Sii 

3 



dVLk'^k^k^ 



An 
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The result for r7(/_free) is 
ijkl _ ■ '^^Ep 



{2nj^-E,)^Tli'' + Tit')+{2nuj-E,-/\YTll'^ , (Bl) 



where the tensor properties are 



rpijkl 

^hh 



l^c,hh / ^c-c- i^XX "^A'X 

m V 20 ' 20 ' ^' 10 ' ^' 



T, 



ijkl 



Ih 



— I Si i^k I H Si kSi / H Si jSj h 

m V 20 ' 20 ' ^' 10 ' ^' 



Tst^ = \\]^{kA,i + kkSj,i) 



and the term involving should not be included if 2hjj < Eg + A. The free particle 
current injection tensor has also been investigated for parabolic bands, but with a simple 
three-band modeL-^^. That model does not have the matrix elements Vc^^ihi and Vd^ih^ and 



thus differs from our result for Tf^^^. 
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